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1.0  INTRODUCTION 

Mesocyclones  are  circulating  regions,  aver¬ 
aging  S  km  in  diameter,  located  within  some  in¬ 
tense  and  well-organized  thunderstorms. 

Although  infrequent,  they  are  an  important  aid 
in  providing  a  storm  warning  signature  detectable 
by  Doppler  radar.  Nearly  all  mesocyclones  are 
associated  with  some  form  of  severe  weather,  and 
about  half  of  them  spawn  tornadoes.  The  typical 
mesocyclone  starts  at  mid-levels  (about  5  km)  in 
a  storm  and  rather  slowly  extends  upward  and 
downward.  Tornadoes  and  other  severe  weather 
manifestations  generally  do  not  occur  until  the 
mesocyclone  base  reaches  the  ground,  a  process 
that  averages  31  minutes  in  duration  (Burgess  et 
al.,  1982).  Consequently,  detection  of  a  meso¬ 
cyclone  early  in  its  development  cycle  will 
yield  a  far  greater  lead  time  for  forecasters 
than  is  presently  available. 

Mesocyclones  have  been  successfully  detected 
for  many  years  through  visual  inspection  of  the 
velocity  field  provided  by  a  single  Doppler 
radar.  Accomplishment  of  this  task  in  real  time, 
however,  demands  the  constant  and  alert  attention 
of  a  meteorologist.  Economical  and  reliable  em¬ 
ployment  of  hvrnian  resources  in  an  operational 
scenario  requires  the  aid  of  an  automated  meso¬ 
cyclone  detection  scheme.  Automation  is  espe¬ 
cially  needed  during  tornado  outbreaks,  when  two 
or  more  mesocyclones  may  occur  simultaneously 
within  range  of  a  radar. 

It  is  proposed  that  this  automatic  proces¬ 
sing  be  performed  by  a  four-dimensional  detection 
algorithm  that  identifies  mesocyclonic  shear  with 
a  resolution-dependent  shear  threshold  applied 
over  a  minimum  velocity  difference.  These  cri¬ 
teria  are  intended  to  enable  the  algorithm  to 
detect  small  mesocyclones,  and  large  ones  at 
great  ranges  (where  Impaired  azimuthal  resolu¬ 
tion  may  reduce  measured  shear)  without  trigger¬ 
ing  false  alarms  due  to  natural  small-scale 
variability  in  velocity.  This  proposed  algorithm 
will  vertically  correlate  mesocyclonic  shear 
features  and  tabulate  their  temporal  persistence. 

The  development  of  this  four-dimensional 
mesocyclone  detection  algorithm  is  a  natural  ex¬ 
tension  of  less  comprehensive  algorithms  pro¬ 
posed  by  Hennington  and  Burgess  (1981)  and  by 
Forsyth  et  al.  (1981).  These  earlier  algorithms 
have  some  complementary  features.  The  former 
detects  mesocyclonic  shear  based  on  a  combination 
of  shear  and  momentum  thresholds,  and  computes 
the  horizontal  extent  of  the  circulation  in 
radial  as  well  as  azimuthal  direction.  The 
latter  scheme  is  based  only  on  shear  between 
velocity  peaks  and  correlatea  this  shear  verti¬ 
cally  and  temporally. 

Momentum  thresholding  was  not  incorporated 
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in  our  mesocyclone  detection  algorithm  because  of 
its  tendency  for  rejection  of  small  but  intense 
vortices,  Buch  as  tornadic  vortex  signatures,  and 
for  acceptance  of  certain  large-scale  patterns 
with  less  than  the  requisite  shear.  We  suggest, 
however,  that  mesocyclones  be  classified  by  the 
trends  of  their  rotational  kinetic  energy,  and 
angular  momentum.  This  may  be  a  means  for  early 
assessment  of  the  destructive  energy  of  the  re¬ 
sultant  tornado.  For  vortices  of  equal  size,  ro¬ 
tational  kinetic  energy  is  proportional  to  the 
square  of  maximum  rotational  velocity,  so  this 
parameter  may  be  especially  sensitive  as  an  indi¬ 
cator  of  the  probability  of  a  tornadic  circula¬ 
tion  penetrating  to  the  ground  and  inflicting 
severe  damage. 

2.0  PROPOSED  ALGORITHM 

The  algorithm  is  capable  of  detecting 
cyclonic  and  anticyclonic  circulations  over  a  wide 
range  of  scales.  It  computes  the  horizonal  and 
vertical  extent,  average  shear,  momentum,  rota¬ 
tional  kinetic  energy  and  temporal  persistence  of 
these  circulations. 

Criteria  for  the  detection  of  mesocyclones 
have  been  set  forth  by  Donaldson  (1970),  and 
later  by  Burgess  (1976) .  Three  basic  require¬ 
ments  are:  1)  significant  azimuthal  shear  must 
exist  between  closed  velocity  contours  of  oppo¬ 
site  sign  (provided  storm  motion  has  been  re¬ 
moved)  .  2)  Shear  pattern  and  closed  isodops  must 

extend  vertically  for  a  height  interval  comparable 
to  the  horizontal  diameter.  3)  The  shear  pat¬ 
tern  must  persist  for  a  time  interval  greater  than 
one  half  of  the  revolution  period  of  the  feature. 

Burgess  et  al.  (1982)  defined  the  evolution¬ 
ary  stages  of  mesocyclones:  1)  The  organization 
stage  is  characterized  by  growth  of  the  circula¬ 
tion  area  upward  and  downward,  from  i.iid-levels. 
Convergence  as  well  as  rotation  may  be  revealed 
in  the  low-level  signatures.  The  organizing 
stage  ends  when  the  mesocyclone 's  base  extends  to 
the  lowest  elevation  scan.  2)  The  mature  stage 
is  characterized  by  a  period  of  maximum  veloc¬ 
ities.  During  this  stage  tornado  potential  is 
greatest.  The  mesocyclone  may  show  some  conver¬ 
gence  at  lower  levels,  pure  core  rotation  at  mid¬ 
dle  levels,  and  some  divergence  at  the  top. 
Generally  the  top  of  the  circulation  is  slightly 
lower  than  the  peak  during  the  organization 
stage.  3)  The  dissipating  stage  begins  when  the 
mesocyclone  summit  drops  rapidly,  and  the  maximum 
velocities  decrease.  At  the  end  of  this  stage 
the  circulation  exists  in  a  shallow  layer  and  is 
very  weak. 


2 . 1  Scale  Considerations 


Although  a  great  majority  of  mesocyclones 
and  tornadoes  are  cyclonic,  a  few  anticyclonic 
tornadoes  have  been  observed.  Therefore,  we 
have  ensured  that  the  algorithm  can  detect  anti¬ 
cyclonic  circulations.  Since  mesocyclones  have 
<ong  been  recognised  as  the  vort icity-producing 
;-.  r-.-e  for  tornad'es  wo  have  included  the  detec¬ 
tion  of  tornado  vortex  signatures  (TVS's)  in  our 
algorithm.  The  algorithm  is  capable  of  detecting 
many  scales  of  meso-circulations .  For  all  prac¬ 
tical  purposes  we  shall  confine  our  analysis  to 
only  those  scales  between  a  small  TVS  and  a  large 
mesocyclone,  that  is,  from  several  hundred  nmtera 
to  20  km  in  horizontal  extent. 

2.2  Algorithm  Processing 

This  algorithm  has  evolved  from  a  two-dimen¬ 
sional  mesecyclonic  shear  algorithm  proposed  by 
Henninqton  and  Burgess  (1981)  (further  developed 
by  Zrnic  et  al.  (1982))  and  a  three  dlsienslonel 
mesocyclone-TVS  algorithm  developed  by  Forsyth 
et  al.  (1981).  We  have  modified  and  merged 
these  algorithms  so  that  they  process  dats  in 
close  to  real  time  and  automatically  archive 
mesocyclone  attributes  for  later  analysis. 

The  algorithm  processing  can  be  broken  into 
five  steps.  These  are: 

U  The  algorithm  searches  for  gradients  in  the 
Doppler  velocity  field.  This  is  accomplished  by 
identifying  either  increasing  or  decreasing  ve¬ 
locities  in  the  a2imuthei  direction. 

2)  When  a  run  of  constantly  increasing  or  de¬ 
creasing  velocities  ends,  the  shear  along  the 
run  is  tested.  If  the  shear  exceeds  a  resolution 
dependent  threshold  a  pattern  vector  is  formed  by 
joining  peaks  it  the  velocity  field  at  constant 
rang*.  The  'mowing  6ix  attributes  of  pattern 
vector  at“  saved:  beginning  and  ending  azimuths 
'■  h$  .  beginning  and  ending  velocities  V.V  ,  range 
to  fhe  pattern  vector,  and  time  of  the  beginning 
az iTj>  h . 

’)  After  each  azimuth  scan  is  completed  the  pat¬ 
tern  vectors  are  grouped  together  into  a  two- 
•  imensiara)  feature.  The  tests  for  association 
of  tv  j  pattern  vectors  are  close  radial  proximity 
(less  that  1  km!  ,  and  azimuthal  overlap.  The 
two-dimensional  features  that  have  a  radial  ex¬ 
tent  of  0.6  ton  (4  range  gates)  or  pore  are  saved 
lor  "rrtical  correlation  with  higher  elevation 
ucanc  .  The  attributes  saved  for  each  two-dlraen- 
lional  feature  are:  the  number  of  pattern  vec- 
to» s  in  the  feature,  minimum  and  maximum  azimuth, 
beginning  and  ending  range,  and  the  minimum  and 
maximum  peik  velocities  and  their  orientation 
-nth  respect,  to  the  radar  beam. 

■‘i  The  vertical  correlation  consists  of  compari¬ 
son  of  the  'cations  of  features  observed  at  suc- 
'■•sr.lve  ole -a*  Ions.  If  the  distance  between  the 
.cotroids  r.f  the  two-dimensional  features  is  less 
than  thr  s  o  <  :  their  diameters  they  are  assumed 
■o  be  correlated  and  their  attributes  are  aver¬ 
age  '  fr*  ■  impart  son  with  features  from  the  next 
pines-  ari.mith  scan.  This  process  continues 
until  fh»  end  of  a  volume  scan  at  which  time  the 
attributes  of  the  three-dimensional  features  are 

’n i 1 v red 

)  V;  recommended  by  Donaldson  (1970)  and  Burgees 
.'19*91  ’•I-  shear  pattern  for  a  mesocyclon  *  should 
persist  r.j  on.-*  hall  the  period  of  the  vortex 
’•ota l  ion  The  algorithm  computes  the  pernls- 
tccc*  of  a  'eature  and  lets  a  flag  when  this 


criterion  has  been  met. 

The  first  test  on  a  three-dimensional  fea¬ 
ture  is  a  check  for  cyclonic  and  anticyclonic 
circulation.  Next  the  azimthal  and  radial  diame¬ 
ters  are  checked,  and  if  they  are  not  within  50 
percent  of  one  another  the  feature  is  laheled  a 
shear  area.  Ctlusiwise,  the  r’i cue-  m  ■ 
pared  to  the  vertical  dep'h  ...  t.  'i.-o  r.e.  ii 
they  are  within  50  percent  of  one  another  the  fea¬ 
ture  is  called  a  mesocyclone.  When  the  shear  for 
a  mesocyclone  exceeds  the  resolution-dependent 
threshold  for  TVS  it  is  labeled  as  a  tornado. 
TVo-dimensional  features  that  are  only  observed  at 
one  elevation  are  called  single-height  features. 

The  results  from  this  analysis  can  inmedi- 
ately  be  written  to  a  screen  for  operator  inter¬ 
pretation,  and  to  a  hard  copy  device.  The  meso- 
cyclone's  average  shear,  momentum,  rotational 
energy,  and  identification  number  are  stored  for 
future  analysis. 

2 . 3  Resolution-Dependent  Shear  Threshold 

In  order  to  compensate  for  the  cross-beam 
degradation  of  resolution  with  range  we  have  ap¬ 
plied  resolution-dependent  thresholds  for  shear 
and  peak  velocity  difference  as  criteria  for  re¬ 
taining  pattern  vectors. 

Hennington  and  Burgess  (1981)  and  Zrnic  et 
al.  (1982)  used  a  combination  of  momentum  and 
shear  thresholds  to  determine  if  a  pattern  vector 
is  saved  for  further  processing.  In  order  to  de¬ 
tect  smaller  circulations  (e.g.,  TVS)  we  have 
eliminated  the  momentum  thresholding  because  it 
implies  a  minimum  length  for  a  pattern  ve-for. 

The  momentum  of  a  pattern  vector  is  estimated 
from  the  product  of  the  measured  velocity  differ¬ 
ence  (V  -V  )  and  the  azimuthal  distance  R (it v-4»e> » 
■hear  ll  the  quotient  of  these  two  quantities. 

The  velocity  pattern  around  a  non-divergent 
mesocyclone  is  generally  modeled  by  a  Rankine 
combined  vortex  (RCV) .  This  means  that  the  core 
of  the  mesocyclone,  defined  as  the  region  between 
velocity  peaks  of  opposite  sign,  is  assumed  to 
rotate  as  a  solid  body.  As  the  mesocyclone  is 
moved  out  in  range  the  peak  velocities  decrease 
significantly  as  more  and  more  of  the  mesocyclone 
is  encompassed  by  the  radar  pulse  volume.  At  the 
same  time  the  diameter  of  the  feature  appears  to 
increase . 

We  have  derived  our  shear  dependence  by  con¬ 
sidering  the  effect  of  increasing  range  of  the 
resolution  of  a  RCV.  Brown  and  lemon  (1976)  pre¬ 
sent  some  enlightening  results  of  a  RCV  model 
developed  by  Zrnic  and  Doviak  (1975).  The  paper 
addresses  the  problem  of  identifying  a  "VS  de¬ 
scribed  by  a  RCV  as  it  is  moved  out  in  range. 

In  this  analysis  we  have  assumed  oata  are 
collected  at  increments  equal  to  the  antenna  i,e  mi- 
width.  The  basic  indicator  of  the  resolution  de¬ 
pendence  of  a  fpature  is  the  ratio  of  the  radar 
beamwldth  (BW)  to  cere  radios  ICR)  .  if  Hi:  «• 
ratio  is  much  less  then  1  the  'ada  -  cm  iim'-e 
the  meBO-circul  at  ion  fait  ly  well.  *».-•  4  *  -*-i  - 

increases  there  is  a  degradation  in  the  observed 
feature.  At  bpamwidth  to  core  ratios  greater 
than  2  the  true  attributes  rf  the  fettmo  -re 
difficult  if  not  impossible  to  resolve. 

r  minimum  aHmuthal  shear  nf  S  -  V-  o  f-» 
mPB-i'-vtlonP  detection  wv.  f-r-i*  j  ' 

Donaldson  (l°7'.l  and  later  corroborated  by 
Butgess  (1976)  from  a  study  of  several  years  of 
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meaocyc lone  data.  If  we  assume  that  an  az¬ 
imuthal  shear  of  5-10"^/s  is  a  practical 
cut-off  between  mesocyc Ionic  and  non-meso- 
cyclonic  storms,  and  that  the  minimum  allowable 
velocity  difference  is  15  m/s,  we  can  see  that 
the  smallest,  lowest  shear  mesocyclone  we  will 
permit  is  3  km  in  diameter.  With  a  beamwidth 
(a)  of  0.8  degrees  we  can  easily  resolve  this 
feature  out  to  a  ranqe  of  110  km  (BW  -  1.5  km) , 
after  which  the  beamwidth  to  core  radius  would 
be  greater  than  1.  At  220  km  our  beamwidth  is 
approximately  3  km  or  twice  the  core  radius  of 
the  RCV.  Brown  and  Lemon's  data  show  that  a 
beamwidth  twice  the  core  radius  would  smooth  the 
peak  velocities  to  approximately  62  percent  of 
their  original  magnitude,  and  the  core  radius 
would  appear  to  be  slightly  larger  than  its  true 
size,  consequently  the  shear  reported  for  the 
feature  would  be  somewhat  less  than  62  percent  of 
the  "true"  shear.  If  the  ratio  of  the  beamwidth 
to  the  core  radius  of  a  tornadic  feature  is 
large,  the  apparent  core  radius  can  be  much 
larger  than  the  actual  size.  For  example,  when 
the  beamwidth  is  3  times  the  core  radius  the 
peak  velocities  detected  by  the  radar  may  be  1.5 
core  radii  from  the  center  of  the  rotation. 

An  illustration  of  the  momentum  and  shear 
thresholds  used  by  Zrnic  et  al.  (1982)  can  be 
seen  in  Fig.  1. 


Fig.  1.  Mesocyclone  detection  thresholds  sug¬ 
gested  by  Hennington  and  Burgess  (1981)  portrayed 
in  terms  of  velocity  difference  and  azimuthal  dis¬ 
tance  across  a  shear  pattern.  H  and  L  are  lines 
of  constant  shear,  and  H  and  L  are  lines  of  con¬ 
stant  angular  momentum.  A  shear  pattern  falling 
into  the  shaded  area  ia  accepted  aa  a  mesocyclone. 
In  region  I,  which  includes  many  high-velocity 
tornado  vortex  signatures,  all  shear  patterns 
would  be  rejected. 


The  shaded  region  of  this  figure  represents  the 
mesocyclone  detection  region.  These  thresholds 
do  not  permit  the  detection  of  TVS's,  as  the  low 
momentum  threshold  precludes  small  but  intense 
vortlcles.  Since  no  objective  criteria  for  the 
detection  of  TVS's  have  been  established  we  Ehall 
use  the  guidelines  adopted  during  the  Joint  Dop¬ 
pler  Operational  Project  (Burgess,  et  al.  1979). 
These  are  the  presence  of  shear  greater  5*10*‘ 
and  a  rotational  velocity  greater  than  30  m/s. 
This  feature  could  have  a  diameter  of  less  than 
1.5  km  and  would  be  in  region  I  of  Fig.  1,  and 
precluded  from  further  processing  with  the  ori¬ 
ginal  momentum  thresholding. 

Fig.  2  is  a  representation  of  our  resolu¬ 
tion-dependent  shear  threshold. 
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Fig.  2.  Mesocyclone  detection  threshold  (in¬ 
cluding  TVS's)  suggested  in  the  present  paper, 
portrayed  in  terms  of  velocity  difference  and 
azimuthal  distance  across  a  shear  pattern.  A 
shear  pattern  observed  with  perfect  azimuthal 
resolution  (Beamwidth/Core  Radius  »  0)  would  be 
accepted  as  a  mesocyclone  if  it  fell  into  region 
I,  With  decreasing  resolution  both  the  shear 
and  velocity  difference  thresholds  would  be  re¬ 
laxed  in  the  direction  of  the  dashed  arrows.  For 
a  Beamwidth/Core  Radius  -  2,  the  mesocyclone  ac¬ 
ceptability  space  is  expanded  to  include  region  I 
and  all  of  region  II. 

The  hatched  region  (II)  of  this  figure  rep¬ 
resents  the  sliding  mesocyclone  detection  region, 
and  the  shaded  region  (I)  represents  the  perfect 
resolution  thresholds.  Using  the  resolution- 
dependent  shear  threshold  one  can  see  that  the 
TVS  discussed  above  would  be  included  in  proces¬ 
sing.  The  shear  threshold  decreases  and  the  min¬ 
imum  velocity  difference  decreases  as  the  BW/CR 
ratio  increases.  The  constant  shear  and  velocity 
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difference  lines  in  Figure  2  repreeent  the  ex¬ 
treme  of  the  thresholds  used  in  the  algorithm. 
Assuming  perfect  resolution  of  a  feature  we  ap¬ 
ply  a  5*i0~3  shear  threshold  and  a  velocity  dif¬ 
ference  threshold  of  15  m/a.  As  the  BW/CR  ratio 
increases  we  gradually  decrease  the  shear  and 
velocity  difference  thresholds  to  2*10"^  and 
?. 1  m/s  at  BW/CR  -  2. 

J.J  MESOCYCLONE  CLASSIFICATION 

Burgess  (1976)  provided  the  earliest  de¬ 
tailed  classif icat ion  of  mesocyclones.  In  the 
first  four  years  of  single-Doppler  data  ac¬ 
quired  by  NSSL,  the  average  difference  in  rota¬ 
tional  velocity  between  tornadic  and  ^on-torna- 
dic  mesocyclones  was  less  than  3ms  .  How¬ 
ever#  the  few  mesocyclones  which  produced  maxi- 
tor  nadoes  had  rotational  velocities  averaging 
more  th<*n  10  m  greater  than  the  non-tornadic 
ones.  These  early  statistics  showed  very  little, 
if  any,  predictive  value  for  tornado  occurrence 
among  mesocyclones*  They  do,  however,  suggest 
that  the  heavily-damaging  max i -tornadoes,  which 
injure  a'*d  kill  ®  disproportionately  large  num¬ 
ber  of  people,  may  be  preceded  by  mesocyclones 
having  unusually  large  rotational  velocities. 

Carrying  this  thought  a  step  further,  it 
seeirs  to  us  that  classification  of  mesocyclones 
according  to  th®ir  rotational  kinetic  energy 
could  be  even  more  useful  than  lotational  veloc¬ 
ity  or  anqclar  momentum.  Rotational  kinetic  en¬ 
ergy  is  proportional  to  the  square  of  velocity, 
so  the  mesocyclones  with  higher  velocities 
woild  be  emphasized  more  than  those  of  larger 
sice.  integration  of  rotational  kinetic  energy 
w*t.t  height  would  provid®  an  energy  index  for 
tad  volume  scan,  and  could  be  an  aid  in 
.ra^ir.q  the  /oluLion  of  a  mesocyclone.  The 
trer.:  zr  tc  i!  r  National  kine.ic  energy  in  a 
me  -.1  cyclcn?  may  enhance  the  ability  of  a  fort- 
ast  <  to  predict  the  touchdown  of  a  maxi-tor- 
nadm . 

4.)  nary  results 

As  e  test  of  the  algorithm  we  have  proc¬ 
essed  several  voljme  scans  of  data  collected 
•»>'» F  (Burmese  et  «1  1579).  Our  teat 

i-..  :=*  arc  fro.r  April  ?0,  1973  starting  at  1010 
CSC,  ju -t  prior  to  the  Piedmont  tornado  touch- 
jnwn .  In  the  first  velum*  scan  the  algorithm 
Identifies  well  davt  loped  mesocyclone.  Out- 
pu  from  the  second  volume  scar,  reveals  that 
ttu  resccyclone  is  moving  to  the  east  and  that 
*t**re  is  a  two-dimensional  TVS  shear  av  eleva- 
.miles  of  5.0  and  6.9  deo-eea,  the  peak 

it  these  levels  being  33  and  -36  m/s. 
wf ,«v.  th  •  .♦r.t9na  ^rops  down  for  the  next  volume 
>  1*  ,  CJT)  we  ha/e  a  wc-dimenstonal  TVS 

I  *ar  wi*h  peak  velocities  of  J5  and  -36  m/a  at 
jnr,  ^  a h»."f  i  v..  jnd  level*  On  the  completion  of 
•i*  t*.  \ !  i  ;r  bun**  scan  the  algorithm  identifies 
*  ional  TVS  . 

in-  ;  i  torradn  touched  down  at  1020 

»:  ,  *\  w!  i:l  t. ioe  the  radar  was  scanning  at 

ejev  t  for.  angles.  The  development  of  thin 
*n-.»*d*  ■-■at  i  ap}  i  at  low  levels,  the  tornado 
.  in  the  j-ound  by  the  time  the  antenna  drop- 
r  *  lo'jrt'pr.  TV: i>  :eat’iie  is  identified 
■  *i  •  :  t  h.v  m  t  w'  more  volume  scans  ar.d 

appe.-**'  weaxen  slightly  although  maintaining 


very  strong  6hear. 

5.0  CONCLUSIONS 

A  resolution-dependent  pattern  recognition 
algorithm  for  detecting  mesocyclones  and  tornado 
vortex  signatures  has  been  developed. 

The  algorithm  identifies  three-dimensional 
cyclonic  and  ant  ley  clou ir  c  ••  • 
putee  their  average  shear,  arg.rai  <..omf;uitrnf  and 
rotational  kinetic  eneigy.  The  time  continuity 
of  the  detected  features  is  confirmed  by  a  persis¬ 
tence  criterion,  and  the  peak  velocities  at  every 
elevation  angle  are  recorded. 

Initial  tests  of  the  algorithm  appear  to  be 
promising.  More  studies  must  be  done  in  this 
area  to  improve  our  understanding  of  mesocyclones 
and  the  mechanisms  involved  as  they  spawn  torna¬ 
does.  The  interactions  between  a  mesocyclone 's 
momentum,  shear,  and  kinetic  energy  as  it  follows 
its  life  cycle  should  yield  some  valuable  statis¬ 
tics  regarding  the  potential  destructive  power  of 
the  circulation* 

Further  refinement  of  the  algorithm  is  like¬ 
ly  to  take  place  as  more  data  is  processed.  In 
particular  the  minimum  feature  size,  the  velocity 
difference  and  shear  thresholds,  and  the  beam- 
width  to  core  radius  ratio  cut-offs  can  only  be 
evaluated  after  extensive  testing.  Probability 
of  detection  and  false  alarm  statistics  need  to 
be  compiled  over  a  wide  range  of  cases  in  order 
to  give  us  better  understanding  of  the  algorithm's 
performance . 
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